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The Similarity Renormalization Group (SRG) is used to soften interactions for ab initio nuclear 
structure calculations by decoupling low- and high-energy Hamiltonian matrix elements. The sub- 
stantial contribution of both initial and SRG-induced three-nucleon forces requires their consistent 
evolution in a three-particle basis space before applying them to larger nuclei. While in principle 
the evolved Hamiltonians are unitarily equivalent, in practice the need for basis truncation intro- 
duces deviations, which must be monitored. Here we present benchmark no-core full configuration 
calculations with SRG-evolved interactions in p-shell nuclei over a wide range of softening. These 
calculations are used to assess convergence properties, extrapolation techniques, and the dependence 
of energies, including four-body contributions, on the SRG resolution scale. 

PACS numbers: 21.30.-x,05.10.Cc,13.75.Cs 



I. INTRODUCTION 

Configuration interaction methods have been used in 
recent years to make increasingly accurate large scale 
ah initio calculations of nuclear structure and reactions 
(e.g., see Refs. [IHS]) Improved algorithms and better 
use of increasing computational resources are critical for 
these successes. However, the reach of these methods 
may also be extended by applying renormalization group 
(RG) transformations to the input Hamiltonian. Renor- 
malization techniques soften the free-space interactions 
by reducing the coupling between high and low momenta, 
leading to improved convergence with the size of the ba- 
sis for a fixed number of interacting nucleons. The Sim- 
ilarity Renormalization Group (SRG) [71 [5] is an attrac- 
tive approach for this purpose due to its relatively simple 
implementation, general flexibility, and the feasibility of 
consistently evolving many-body operators [21 HOj . 

Previous studies of the SRG in nuclear physics estab- 
lished its usefulness for few-body systems by demonstrat- 
ing improved convergence with two-nucleon (NN) inter- 
actions alone [TTHTB] . In Ref. [H], a detailed study of 
SRG convergence with NN forces in the p-shell was made. 
The present work extends this study to include initial and 
induced three-nucleon (NNN) forces. We build upon the 
technology to evolve NNN forces introduced in Ref. [T5] . 
which was applied in Ref. [T^ to "^He and ^Li to ex- 
plore the effects of full two-plus-three-body evolved in- 
teractions in light nuclei. Roth ct al. have subsequently 
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used the Importance Truncated No-Core Shell Model 
(IT-NCSM) with the SRG [H] to significantly push the 
limits in A and model space size [IHl HI] . We focus here 
on a wider set of p-shell nuclei and on wider ranges of 
softening that, with our extrapolation methods, yield 
nearly converged results without implementing impor- 
tance truncation. 

The SRG flow equations generate a continuous se- 
ries of Hamiltonians that is analogous to the running 
of the strong coupling constant in quantum chromody- 
namics, having both scale (or resolution) and scheme de- 
pendence [lOj . The scheme dependence arises from the 
choice of initial nuclear Hamiltonian and the choice of the 
operator generating the flow (see below) . While the SRG 
offers a useful framework for future comparisons of such 
choices for both NN and NNN interactions and exploring 
the flow to universal forms [HI [501, in this work, we restrict 
our attention to just one choice. In particular, we use the 
chiral effective field theory (EFT) potential at N'^LO with 
500 MeV cutoff from Ref. [21] together with an NNN po- 
tential at N^LO 122] in the local form of Ref. [23,. This 
is also the Hamiltonian used in Refs. [THl [1^1 [T51 [21] . 

The scale dependence arising from the RG flow is man- 
ifested as a decreasing decoupling scale that marks the 
energy difference for which matrix elements between off- 
diagonal energy states become highly suppressed. For- 
mally, all of the evolved Hamiltonians have equivalent 
physics content to that of the initial Hamiltonian, so it 
would seem to be advantageous to evolve to very low 
scales to optimize the convergence of many-body calcu- 
lations. But in practice, the initial and running Hamilto- 
nians are expanded in a finite basis (harmonic oscillators 
here) and many-body forces are truncated at some level. 
Therefore, it is necessary to monitor and characterize the 
evolution of many-body forces and the residual running 
of calculated observables, which can vary with the size 
of the nuclear system. In this paper, we present bench- 
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mark calculations using SRG-evolved interactions in p- 
shell nuclei. We use them to explore the characteristics 
and practical limits of SRG evolution for these systems 
by assessing convergence properties, extrapolation tech- 
niques, and the stability of predicted observables. 

In Section |llj we briefly review the formalism used in 
this study and summarize some observations from previ- 
ous work. A more complete discussion is provided in |16| . 
The convergence properties of the evolved NN and NNN 
potentials in various nuclei are explored in Section |III| 
and in Section IIVI we examine the evolution with A and 
the patterns of induced many-body forces. Section |V] 
summarizes our conclusions and provides an outlook for 
future calculations. 



TABLE I. Definitions and values of basis-size truncations used 
on the initial Hamiltonian and on the model space for the 
CI calculations. The value of A'^ASmax is uniformly reduced 
for three-body partial-wave JttT channels with higher J to a 
minimum of 20. 

AA2max Maximum sum of two-body oscillator 300 
quanta for initial two-body matrix element 
evolution 

AA3max Maximum sum of three-body oscillator 40 
quanta for initial three-body matrix ele- 
ment evolution 

A^max Maximum sum of j4-body oscillator quanta 2-8 
above the minimum for the A-body system 



II. BACKGROUND 

A. SRG evolution 

The SRG for low-energy nuclear physics generates a 
continuous series of Hamiltonians Hx from an initial free- 
space Hamiltonian H\^oo by unitary transformations, 

Hx = UxHx=ooUl , (1) 

which is carried out by solving a set of flow equations, 

(2) 

With an appropriate choice of the hermitian operator 
G, the Hamiltonian is evolved to band-diagonal form 
with respect to energy [HI [25]. In most nuclear appli- 
cations to date, the relative kinetic energy Ti-d has been 
used because of the favorable convergence properties of 
the evolved Hamiltonians and for its convenience in con- 
structing basis expansions [H [10] . We use it here exclu- 
sively but note that other choices may be advantageous 
in reducing the growth of many-body forces. This is be- 
ing investigated separately [261 EZ] • The flow parameter 
A keeps track of the sequence of Hamiltonians (s or a, 
with s = a = 1/A^, are also used elsewhere [TTl [T^ fT5]). 
For G = Tici, A has dimensions of momentum and runs 
from oo toward zero with increasing softening. 

Evolution is performed in the Jacobi-coordinate har- 
monic oscillator (HO) basis used for the No-Core Shell 
Model (NCSM) [2[H]. This is a translationally invariant, 
anti-symmetric basis for each A-body sector, in which a 
complete set of states in the model space is defined by 
the maximum excitation of N^^i^-^hfl above the minimum 
energy configuration, where fl is the harmonic oscilla- 
tor parameter. This basis is variational in A^max; that 
is, the energy converges asymptotically from above as 
more basis states are included. The SRG preserves this 
variational characteristic through smooth unitary evolu- 
tion, in contrast to Okubo-Lee-Suzuki based renormal- 
izations [21IS], which are unitary transformations specific 
to the model space. 



We start by evolving H\ using Eq. Q in the A = 2 
subsystem, completely fixing the evolved two-body ma- 
trix elements. Next, by evolving H\ in the A — 3 subsys- 
tem we determine the combined two-plus-three-body ma- 
trix elements. We isolate the three-body matrix elements 
by subtracting the evolved two-body elements within the 
A — 3 basis [35]. Having obtained the separate NN and 
NNN matrix elements we can apply them as inputs to 
any ab initio nuclear structure problem. We are also free 
to include an initial three-nucleon force in the starting 
Hamiltonian without changing the procedure. 

While any initial interaction can be used as input to 
the SRG evolution, here we use the chiral EFT NN po- 
tential from the 500 MeV N^LO interaction of Ref. [2l] 
exclusively. As an initial NNN potential, we use the chi- 
ral N^LO potential [21] in the local form of Ref. [35] . 
The low-energy constants cd — —0.2 and ce — —0.205 
are the result of a fit to the average of triton and '^He 
binding energies and to triton beta decay as described in 
Ref. [33]. 

Hamiltonians obtained via free-space SRG evolution 
are independent of the basis choice if the basis is suffi- 
ciently complete. That is, a Hamiltonian evolved to a 
given A reproduces the results of a Hamiltonian evolved 
to the same A in a different basis. But in practice there 
are truncations, both in A-body forces and in basis size, 
that are relevant to controlling the quality and consis- 
tency of SRG-evolved interactions. 

In the present work, induced four-body (and higher) 
forces are not included, so calculations for A> 4 will be 
only approximately unitary. The many-body interaction 
matrix elements induced by the evolution appear in a de- 
creasing hierarchy in few-body nuclei |15l I16j . One of our 
goals is to determine if that hierarchy is maintained for 
p-shell nuclei or if the induced many-body contributions 
become unnaturally large for certain systems and/or val- 
ues of A. 

Because of computational constraints, we are forced to 
apply separate truncations to the A = 2 and ^ = 3 sec- 
tors of the initial Hamiltonian, which we denote A^A2max 
and iVA3max, respectively (see Table |l|. These cutoffs in 
the basis size must be large enough to fully accommodate 
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the ultraviolet (UV) contributions (or high-momentum 
components) from the initial NN plus NNN Hamilto- 
nian. The ultraviolet cutoff in an oscillator basis scales 
like \/NMl, where N represents the maximum number of 
single-particle oscillator quanta in the basis, so there will 
be an Ml below which the initial and, therefore, evolved 
Hamiltonian projections onto the oscillator basis are in- 
complete. When we use such an Ml that is too low, 
we are effectively working with a different Hamiltonian. 
As a consequence, the calculations of observables in the 
many-body basis with the too-low Ml will not converge 
(or extrapolate) to the same results found at larger Ml. 
This is not a problem for the NN interaction, for which 
-^A2max IS Sufficiently large for the chiral EFT Hamilto- 
nian for all Ml considered, but becomes a factor for the 
NNN force, as illustrated below. 

Hamiltonians are derived and evolved in the Jacobi 
basis for A = 2 and 3 and then translated to a Slater de- 
terminant basis for full configuration interaction (CI) cal- 
culations of larger systems. The particular CI procedure 
used here, including the extrapolation to infinite basis 
size and associated uncertainty estimates, is referred to 
as no-core full configuration or NCFC [T]. Other CI cal- 
culations in the literature using SRG-evolved interactions 
and extrapolation are called NCSM (e.g., Refs. P^fTB] ). 
While the original NCSM featured a finite matrix trun- 
cation and an effective Hamiltonian renormalized to that 
finite space [5], these SRG-based NCSM and NCFC pro- 
cedures are equivalent except for variations in the extrap- 
olation and uncertainty quantification procedures. 

In these CI calculations, the size of the largest feasible 
model space is highly constrained by the total number of 
two- and three-body matrix elements in the full space. 
Fortunately, the MFDn code [3DH32] that carries out the 
Lanczos matrix diagonalization algorithm is highly opti- 
mized for parallel computing. The calculations were per- 
formed on the Intel Xeon cluster Sierra at LLNL, using 
up to about 15 TB of memory across 7,200 cores on the 
Cray XE6 Hopper at NERSC, using up to about 100 TB 
of memory across 76,320 cores; and on the Cray XK6 
Jaguar at ORNL, using over 500 TB of memory across 
261,120 cores. MFDn has been demonstrated to scale 
well on these platforms for these types of runs |33| . 



B. Guide to the calculations 

To set the stage for our examination of p-shell nuclei, 
we show in Fig. [T] the ground-state energy of "^He as a 
function of A [TSl [T5]. Here, and throughout the pa- 
per, we will compare three types of calculations, which 
are summarized in Table |ll) The NN-only calculations 
include two-body matrix elements that are phase-shift 
equivalent to the initial NN interaction. When this part 
of the Hamiltonian is used alone to study A > 2 sys- 
tems, the results are not unitarily equivalent at differ- 
ent A because the SRG has reorganized the degrees of 
freedom to reduce coupling of high- and low-energies. 



TABLE II. Guide to the calculations. 

NN-only No initial NNN interaction and do 
not keep NNN-induced interaction. 

NN + NNN-induced No initial NNN interaction but keep 
the SRG-induced NNN interaction 
arising from the NN interaction 
alone. 

NN + NNN Include an initial NNN interaction 
and keep the SRG-induced NNN in- 
teraction arising from the combina- 
tion of NN and NNN interactions. 




X [fm 

FIG. 1. (color online) Ground-state energy of ''He as a func- 
tion of A for the three calculations in Table HIl The results at 
all A are converged at the 2 keV level or better. The dotted 
line is the unevolved result. 



Formally, in the process of maintaining unitary equiv- 
alence in all sectors the SRG induces new contributions 
to many-body matrix elements, but in the NN-only re- 
sult these induced interactions are omitted. We get a 
characteristic pattern (see the NN-only results in Fig. [T]) 
where the converged result varies with evolution, starting 
at some (underbound) level for the initial Hamiltonian, 
then falling and rising again with subsequent evolution 
(decreasing A). 

In the NN+NNN-induced calculations, the Hamilto- 
nian includes the evolved NN matrix elements as well as 
all three-body matrix elements induced by the SRG start- 
ing from only an NN interaction (i.e., no initial NNN). In 
an ^ = 3 system this will be unitarily equivalent to the 
initial NN-only Hamiltonian, so the energy spectrum will 
be the same, up to numerical truncation errors (e.g., be- 
cause of an insufficient A^ASmax)- Finally, the NN-I-NNN 
calculations include an initial three-body interaction as 
well as the induced three-body matrix elements that now 
arise from the combined evolution of the NN-f NNN in- 
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teractions. For A = 3, it should be unitarily equivalent 
to the initial NN+NNN Hamiltonian. For A > A, there 
will be induced four-body (and higher-body) interactions 
that are not included in any of the present calculations. 
Their omission causes differences in the predicted energy 
spectra as a function of A. 

The computed ground-state energies for ^He in Fig. [l] 
are well converged at all A, so the interpretation is clear. 
For the NN-only calculations, deviations from unitary 
equivalence are evident just below A = 10 fm""'^, where 
binding is increasing by a maximum of about three MeV 
(10% of the total binding energy), peaking just below A — 
2fm"^ and then decreasing rapidly and overshooting the 
original ground-state energy by A = Ifm^^. A similar 
pattern for NN-only was shown in Ref. |14j for several 
p-shell nuclei. 

The NN-f NNN-induced calculation shows a dramatic 
reduction in the variation of the energy for A > 1.5 fm~^, 
with only a small decrease in the binding energy peaking 
near A = 2fm~^. The deviations near A = 1 fm~^, which 
imply net induced four-body interaction contributions, 
are only about 300 keV, or still an order of magnitude 
reduced from the largest NN-only variations. The same 
pattern for the (implied) induced four-body interaction 
is seen when an initial NNN interaction is included, with 
just a slight change in the pattern at A = Ifm""'^. When 
we compare to the larger nuclei in the present work, we 
will not be able to examine the full range of A used in 
Fig. [T] because convergence is only sufficient for reliable 
extrapolation with small errors for A up to about 2fm^^. 



C. Extrapolation methods 

For well-evolved Hamiltonians in lighter nuclei (e.g., 
or '^He at A < 1.5 fm~^), our predictions for ground- 
state energies are well converged at computationally ac- 
cessible values of iVmax- However, for larger nuclei and 
greater A values we will need to extrapolate calculated 
energies to iVmax = oo. To do so, we primarily use em- 
pirical extrapolation procedures based on those described 
in Refs. [TJ [TH [TB] but also compare to a new procedure 
from Ref. 

The empirical model used for ground-state energies is 



Ar 



(3) 



where Aa and ba are {Ml. dependent) constants, Nai are 
the A^max values, and a labels the HQ value. The goal 
is to determine the common parameter E^^, which is 
the estimate for the ground-state energy extrapolated to 
A^max = oo. This can be cast as a one-dimcnsional con- 
strained minimization problem with the function 

g{Eoo) = ^(log(£;„ - E^) - a„ - b^N^if/al^ , (4) 



where the {aa} and {6q} are determined directly within 
the function g by invoking a constrained linear least- 
squares minimization routine. The constraint is the 



bound Eoo ^ min{{Eai}), where E^o < and "min" 
means "most negative" . (One can also allow for weights 
depending on A^ax and/or hil.) 

In the present investigation, we have applied this ex- 
trapolation model for individual values of hfl, determin- 
ing error estimates as in "Extrapolation B" from Ref. [T] , 
but also including several values of HQ in a constrained fit 
over a range where they are considered reliable. We em- 
phasize that while this model has been generally success- 
ful when applied in NCFC calculations with SRG-evolved 
interactions, its validation is empirical rather than theo- 
retical. 

An alternative EFT-motivated approach to extrapo- 
lation is based on explicitly considering the ultraviolet 
(UV) and infrared (IR) cutoffs imposed by a truncated 
harmonic oscillator basis [Ml [3S] . This has led to a the- 
oretically motivated IR correction formula and an em- 
pirical UV correction formula |34| in which the basic 
extrapolation variables are the effective hard-wall size 
L and the analogous cut-off in momentum, Aijv- In 
terms of the oscillator length b = yJh/{mQ), rough esti- 
mates of these variables are Kuv ~ ^/2{N + 3/2)h/b and 
L « ^/2iN + 3/2)b, where A^ = A^,„ax + 1 for p-shell nu- 
clei [311 A formula combining both corrections (i.e., 
they are treated independently) takes the form [34 



E{Auv,L)^E, 



(5) 



Note that this formula contains exponentials with argu- 
ments proportional to both A^ (from Af jy) and \/N (from 
L), in contrast to Eq. ([3|. 

Following Ref. [M], we apply Eq. ^ with Erx>, Bq, 
Bi, B2, and fcoo treated as fit parameters that are deter- 
mined from a simultaneous optimization to data at all 
Ml^ including the intermediate region where both IR and 
UV corrections are significant. It may be advantageous 
in general to isolate the IR or UV corrections by using 
only large Kl or small hVl results, respectively. How- 
ever, most of the present calculations were made with Ml 
values close to the energy minimum, which means com- 
parable UV and IR contributions [34]. (The exception 
is for very low A, where UV convergence is reached for 
all Ml considered.) We also exploit a recent observation 
that the expressions for L and Kuv give more accurate 
energy corrections if we take N N + 2, which is par- 
ticularly effective when A^ax is small 36J. Thus we will 
use A^ = (A"m ax + 1) + 2 = A^max + 3 for the calculations in 
Section HI C Equation ^ has been successfully applied 
to NN-only calculations from Ref. [14], but here we test 
it for the first time with three-body forces included. 



III. CONVERGENCE 

For fixed Amax, both the UV and IR momentum cut- 
offs scale with y/Ml, which means that there is a trade- 
off: increasing Ml increases the ability to accommodate 
high-momentum components while decreasing the ability 
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to accomodate long-distance physics. The result is a fa- 
miliar variational minimum with respect to Ml. For NN- 
only calculations, it was observed in Ref. [13] that with 
decreasing SRG A at fixed A^max, the minimum systemat- 
ically shifts to lower hfl and convergence becomes much 
more rapid. Here we examine if these observations are 
modified by the presence of a three-nucleon force (3NF) 
in NN-KNNN-induced and NN-^NNN calculations. 



A. Size of three-body evolution basis 

As noted earlier, the size A^ASmax of the A = 3 ba- 
sis we use to evolve the Hamiltonian before embedding 
in larger systems is limited by computational constraints. 
In Fig.[2j the impact on the calculated ground-state of the 
triton is shown for the unevolved interaction using differ- 
ent values of A'Asmax with hfl ranging from 10 to 24MeV. 
At each HQ, the signal that iVAamax is sufficiently large is 
convergence of the ground-state energy, which is evident 
for hfl > 18. In contrast, the systematic underbinding 
at lower hfl values with the iVASmax from Table |l] will 
be preserved when the Hamiltonian is evolved; in effect 
a different initial Hamiltonian will be used. The spread 
of points at fixed HQ and particularly the deviations at 
TVASmax = 40 from the fully converged energy imply that 
low oscillator parameters (i.e., below hil = 18MeV), will 
be unreliable for energy calculations in larger nuclei (be- 
cause we cannot predict the degree of underbinding, as 
shown in Fig. |4] discussed below). In contrast, the size 
of the A — 2 basis used here is sufficient for convergence 
within 1 kcV in the full range of hft considered. 
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FIG. 2. (color online) Ground-state energy of the triton 
for the unevolved chiral EFT Hamiltonian in different three- 
body basis sizes (A'^ASmax) with a large, fixed two-body basis 
{Na2 max — 300). 

The implications of a too-small iVASmax sA, smaller 
for SRG evolution are illustrated in Fig. [3j which shows 
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many-body basis size of A^max ~ 8 for Hamiltonians evolved 
to A = 1.5 fm^^ in different three-body basis sizes (A'ASmax) 
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fixed Amax = 8 of *Be ground-state energies from Fig. [s] at 
several values of HQ. 



fixed A'lnax = 8 (and fixed A^A2inax — 300) calculations of 
the ^Be ground state of the same initial NNN interaction 
evolved in different three-body basis sizes, for a range 
of hfl. All results are for A = 1.5 fm~^. It is evident 
that hil. = 12 MeV is far from converged even at the 
largest value of A^ASmax available. Note that we do not 
expect the different hfl calculations in this figure to have 
the same energy as A'ASmax — ^ oo, because Af^ax = 8 is 
still unconverged in the many-body system. However, for 
each individual Ml we need convergence for the largest 
A^ASmax, as obscrved for hfl — 20 MeV. We also note 
that an exponential model for the convergence in A^ASmax 
does not work, as shown in Fig. |4] For hft = 20 MeV we 
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observe good convergence at iVASmax = 40. However, the 
quality of the exponential fits in Fig. |4] deteriorates as 
fi.r2 decreases below 20MeV. We conclude that simple 
exponential extrapolation in iVASmax is not an option for 
hft — 16MeV and below. In particular, while hH. — 
16MeV is close to converged for ^Be, this may be less 
true for larger nuclei, so we will (perhaps conservatively) 
only consider fiO > 18 MeV to be reliable in the following. 



B. Convergence with model space size N„: 



In Figs. [5] through 13 we show detailed results for the 
ground-state energy for each of the calculations of Ta- 
ble |IT] as a function of hft for three representative nuclei 
(^Li, ^"B, and ^^C) and basis sizes A^max from 2 to 8 
(and 10 for some NN-only cases). Panels from left to 
right show results for decreasing SRG A. In all panels 
the symbols connected by dashed lines denote the ab ini- 
tio calculated points at a single iVmax, while the solid 
lines are exponential extrapolations at each individual 
hft based on "Extrapolation B" from Ref. [1 . Extrapo- 
lations are given for A^max — 2-6 and A^max = 4-8, with 
error bars in the latter determined by the difference from 
the central prediction of the former. For ^Li and ^"B 
NN-only calculations, there are also extrapolations for 

4-8; see 
= 10 



6-10 with error bars based on N„ 



Fig. 5 of Ref. 



for C for results up to N„ 



with the same NN-only interactions, but without the 
Coulomb interaction. These figures display the system- 
atics of NCFC convergence of SRG-evolved interactions 
with and without three-body forces. 

Here, we make various summary observations based on 
these figures (with further discussion of extrapolations in 
the next section): 

• For all calculations, the rate of convergence with 
basis size iV,„ax is greatly accelerated with decreas- 
ing A for fin near the variational minimum. The 
addition of NNN interactions does not affect the 
rate substantially. The dependence on hSl at low A 
is very flat for higher A^max- 

• In almost all cases, the location of the variational 
minimum in Ml for a given A^max shifts to smaller 
hH. as A decreases. This is expected because the 
minimum is where the ultraviolet (UV) and in- 
frared (IR) corrections are roughly equal [31] . (See 
Fig. [17] and the accompanying discussion below.) 
The contribution from high momentum compo- 
nents of the interactions to low-energy states de- 
creases as A decreases, so lower A Hamiltonians 
are less sensitive to the UV reach in momentum 
of a truncated HO basis. In other words, the more 
evolved Hamiltonian is less sensitive to the UV cut- 
off intrinsic to the HO basis and thus the inter- 
play between UV cutoff and IR distortion is shifted. 
This results in a systematic shift of the optimal hfl 
to lower values as the Hamiltonian is evolved. 



• The location of the minima for NN-|-NNN-induced 
calculations are systematically lower in hfl than for 
the corresponding NN-I-NNN interaction calcula- 
tions. This is consistent with the induced NNN 
interaction being softer than the initial NNN inter- 
action. The significance of high momentum com- 
ponents in the initial NNN interaction were already 
evident above where we saw that they are not as 
well converged in the A — 3 space used for evo- 
lution and require larger Ml. Furthermore, with 
NN-|-NNN-induced interactions, the nuclei are sys- 
tematically less bound than with the NN-I-NNN in- 
teraction, so the wavefunctions will have longer- 
range exponential tails; as a consequence, the cor- 
responding wavefunction is better represented with 
a lower hil for NN-|-NNN-induced calculations than 
for NN-I-NNN calculations. 

• Results for A = 1.5 fm^^ and 1.0 fm""'^ are suffi- 
ciently converged that extrapolations in the non- 
gray regions (gray shadowing signifies the unreli- 
able region of hil discussed above) are the same for 
different hfl within the (small) error bars. Near 
the minima in this A range extrapolation is largely 
superfluous. In cases including initial or induced 
NNN interactions, the minimum is sometimes in 
the gray region. However, in such cases the conver- 
gence for larger hfl is well advanced and the indi- 
vidual Ml, extrapolations are consistent with each 
other. 

• Extrapolated results for A = 2.5 fm~^ (and to some 
degree for A = 2.0 fm~^) appear to depend sys- 
tematically on hil. Although it has generally been 
considered most reliable to extrapolate using ener- 
gies for hfl near the variational minimum, selected 
results in A'^max = 10 spaces (e.g., for -^°B NN-only) 
suggest that the extrapolations from 4,6,8 are over- 
bound near the minimum for larger A (and 2,4,6 
even more so) , so that estimates from larger fifi are 
more robust. This may be related to the fact that 
the location of the variational minimum is shifting 
to larger values of fifi as A^max increases. In general, 
our A > 2.2 fm~^ extrapolations for these nuclei 
have error bars too large to allow useful quantita- 
tive conclusions about A dependence. 



• In Figs. 



11 



12 



and 



13 



the A = 2 fm results 



for ^^C at hil — 20 MeV include a horizontal line 
marking the best extrapolated value (with error 
bar) from the IT-NCSM calculations of Roth et 
al [IB]. The A'^max < 6 values are relatively con- 
sistent with the extrapolated IT-NCSM results, but 
beginning with the extrapolated NCFC results that 
include the A'^ax = 8 points, we appear to predict 
somewhat more binding than the IT-NCSM extrap- 
olated result, though our uncertainties are large 
enough that a definitive conclusion is elusive. Note, 
however, that the prescription for the A^ASmax trun- 
cation is slightly different in Ref. [T5] than what 
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we use here. Both in Ref. 



and in our calcula- 



tions use A^ASmax = 40 for the leading, J = 5, 3NF 
contributions, but we use different cutoffs for the 
higher- J terms. A detailed analysis of the effect of 
importance truncation for this particular case is in 
progress |37^. 

• It is evident by comparing extrapolated values for 
NN-only to the NN+NNN calculations that the A 
dependence of the extrapolated energies is signifi- 
cantly reduced by including 3NFs. This is summa- 



rized in the figures discussed in Section IV 



Overall, the convergence patterns with iV,„ax and Ml pre- 
viously observed as a consequence of NN-only SRG evo- 
lution [M] are still present when 3NF contributions are 
included. 



C. Comparison of extrapolation methods 

The extrapolations in Figs . [s] through [T3| are a series of 
individual extrapolations using Eq. ([3]), each for a fixed 
value of hfl. As noted earlier, the plotted error bars 
for each hfl are determined by comparison of each indi- 
vidual A'niax = 4-8 extrapolation to that obtained from 
A'max = 2-6. The resulting set of predictions can then be 
analyzed to obtain a predicted energy and overall error 
bar [1]. We consider two criteria for determining a best 
"Extrapolation B" result for Eao and the corresponding 
error bar: 

• the result at the hfl value for which the amount 
of extrapolation is minimal (i.e. the point where 
E{Nmi,x) - Eoo is minimal); 

• the result at the hfl value for which the numerical 
error estimate is minimal. 

In both cases, we restrict the best "Extrapolation B" to 
hfl values at or above the variational minimum at the 
highest iVijiax employed in the extrapolation. For the 
error estimates we use the average of the error bars in a 
region of 8 MeV around this best hfl value. This initial 
error estimate is enlarged as necessary in order to get 
consistent results, such that the central values are within 
the error estimate in the entire 8 MeV range. 

Note that below A = 1.5 fm"^ with 3NFs, we cannot 
apply either of these criteria, because we are only using 
results for hfl > 18 MeV, and the variational minimum is 
typically at or below hil = 18 MeV for the lowest A val- 
ues; however, for these values of A the results are quite 
close to convergence, and we base our error estimate on 
the results for 18 < < 24 MeV. Above A = 1.8fm~\ 
the two criteria give (slightly) different results, but gen- 
erally with overlapping extrapolation error estimates. 

The convergence pattern and extrapolations of the NN- 
only data up through iVmax = 10 (see Figs. [5] and [s] as 
well as Fig. 5 of Ref. |33]), suggest that for A — 2.5 fm^^ 



the hfl value that minimizes the numerical error esti- 
mate, hn ~ 36 MeV, is more reliable for the extrapo- 
lation (at least at this A value) than the hfi value that 
minimizes i?(iV,nax) — E^o- (With the nonlocal NN in- 
teraction JISP16, for which this extrapolation has been 
used extensively [U [3l |38j these two criteria rarely lead 
to significant differences.) 

An alternative approach is to use a constrained opti- 
mization that uses all (or a specified subset) of hft data, 
requiring the same extrapolated energy i5oo iu Eq. ([3| for 
every hil, such as "Extrapolation A" of Ref. [T|. Here, 
we use basically the same procedure for the constrained 
optimization, using four subsequent iVmax values and five 
hn values for each constrained fit, with the same weights 
and error estimate as in Ref. |T] . (For the current calcula- 
tions five hfl values span an 8 MeV range of hft, whereas 
in Ref. Ij it spans a 10 MeV range.) Again, with 3NFs 
this procedure cannot be applied below A = 1.5 fm~^, but 
the results are close to convergence for these cases, and 
therefore less sensitive to the details of the extrapolation. 
For 1.5 fm~^ > A > 2fm~^ the results from "Extrapola- 
tion A" are consistent with those from "Extrapolation 
B", but for A > 2fm-^ the procedure of Ref. [T] for "Ex- 
trapolation A" , modified as described above, leads to re- 
sults that show a systematic deficiency at A'^^ax = 10. 



• x = 


= 2.5, ha = 26 


♦ X = 


= 2.5, hO. = 36 


dashed: 


single 2-6 


dotted: 


single 4-8 



> 




FIG. 14. (color online) Ground-state energy extrapolations of 
'^Li as a function of A^max with an N'^LO NN interaction |21j 
evolved to A = 2.5 fm~^. The symbols are the calculated 
points. The curves show single extrapolations using Eq. (|3f 
with iVmax = 2-6 (dashed), 4-8 (dotted) and 6-10 (solid) at 
(blue) ftfi = 26 MeV which minimize the amount of extrapola- 
tion and at (red) HQ — 36 MeV which minimize the numerical 
error estimate. The horizontal dotted and solid lines, with 
the band indicating the associated error bars, are the result 
from a constrained fit following the procedure of Ref. il; for 
five hn values from 22 to 30 MeV. 

Figure [M] shows a comparison of the single and the 
constrained extrapolation schemes for A = 2.5 fm""'^. In 
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this example we use '^Li with NN-only calculations up to 
-^max = 10. The variational minimum is at 24MeV, so 
we use five hil values from 22MeV to 30MeV, and four 
-^max values for each h^l for the each of the constrained 
fits. Figure 14 shows clearly that this procedure with 
the Njnax = 2 to 8 results leads to an overestimate of 
the binding energy and an underestimate of the extrap- 
olation uncertainty. In Fig. [T4j we can also see why this 
procedure leads to erroneous results: near the variational 
minimum the convergence is not a simple exponential for 
this value of A, as is evident from the single hH. fits (blue 
curves) . 
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FIG. 15. (color online) Ground-state energy of ^Li for 
NN-I-NNN evolved Hamiltonians at A = 2.0 fm~^, plus extrap- 
olations based on "Extrapolation B" from Ref. [T] (solid line 
and points with error bars) and on the combined IR/UV cor- 
rection formula Eq. Q that yields Eoo (dashed line) and indi- 
vidual corrections for each hO, and A^max combination (points 
near the dashed line) based on the single set of best-fit pa- 
rameters. 

The alternative EFT-motivated approach described in 
Sec. |II C| is complementary to the schemes used in Figs. [5]- 
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A sample application of Eq. ^ , which includes both 
UV and IR corrections to the energy at each iVmax and 
h^, is shown in Fig. [15] for the NN-I-NNN calculation of 
the ^Li ground state. In this example, the fifteen points 
with hil > 18 MeV and A'max from 4 to 8 are inputs to 
a simultaneous fit of the five parameters of E{Auv,L) 
with AjjY and L given as functions of Ml and A'max 



(see Sec. II C for the formulas, recalling that we will use 
N — A^max + 3 for p-shell nuclei). The result for Eoo 
(dashed line) is consistent with the individual "Extrap- 
olation B" results (from Fig. [?]) and the overall result is 
within the error bars of that scheme. An error analy- 
sis procedure for the IR/UV correction model is not yet 
available, which limits its utility for the present analy- 
sis, but the goodness-of-fit can be assessed by the (very 
small) spread of corrected points about the dashed line. 
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FIG. 16. (color online) Ground-state energy of for the 
NN-I-NNN evolved Hamiltonians at A = 2.0 fm"'^, plus extrap- 
olations based on "Extrapolation B" from Ref. [T] (solid line 
and points with error bars) and on the combined IR/UV cor- 
rection formula Eq. Q that yields E^c (dashed line) and indi- 
vidual corrections for each hi} and Vmax combination (points 
near the dashed line) based on the single set of best-fit pa- 
rameters. 



Corrected points represent the use of the best fit parame- 
ters with Eq. ^ (except Eaa) to extend each finite basis 
result to infinite iVmax at fixed fifl. Note also the predic- 
tions for the hil — 16 MeV points, although those data 
points were not included in the fit. 

A second extrapolation based on the IR/UV correc- 
tion model is shown in Fig. 



16 where the spread of cor- 



rections indicates a still good but less-than-ideal fit for 
the NN-^NNN calculation of the ^^C ground state. Note 
that the largest deviations of the corrected results from 
the fit Eoo occur for two of the three highest TVmax points. 
The implication is that the true Eoo should be slightly 
more negative, which is also the conclusion from compar- 
ing with the "Extrapolation B" analysis. Note, however, 
that although we used the same total number of points for 
this IR/UV extrapolation, we have only three A^max = 8 
points (see Fig. 13 compared to five of the fifteen points 



that we used for Fig. 15 



In general, a good fit requires the UV and IR func- 
tional forms to be adequate models for smaller A^max 
values (with increasing A^max there is decreasing sensi- 
tivity while the computational cost is increasing dramat- 
ically). Much remains to be explored for heavier nuclei 
but detailed investigation of two-particle models and the 
deuteron suggest that the forms in Eq. ^ can be im- 
proved and that the predictions can be sensitive to op- 
timizing the choice of expressions for L and Ayy |36j 
(e.g., using A^^ax + 3/2 -I- 2 rather than A^^ax + 3/2 for 
two-body systems). We are also not yet able to take ad- 
vantage of the theoretical prediction that koo should be 
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related to the nucleon separation energy and the empiri- 
cal observation that Bi is found to be numerically close 
to A. 

The two examples considered so far are for A = 2 fm~^. 
For A < 1.5fm~\ the new extrapolation method using 
data up to A^max — 8 gives predictions for the iV,„ax — oo 
energies consistent with the other extrapolation schemes. 
As already noted, there are some systematic differences 
for A = 2fm^^, but they are within the "Extrapolation 
B" uncertainties. For larger A these differences grow, but 
it is not possible at present to determine which approach 
is superior. 
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FIG. 17. (color online) Ground-state energy of ^Li for the 
NN-I-NNN evolved Hamiltonians at A = 2.0fm"\ with IR 
(vertical dashed) and UV (vertical dotted) corrections from 
Eq. ([5| that add to predicted Eoc values (points near the 
horizontal dashed line, which is the global Eao)- 



It is instructive after making a global fit to decompose 
each correction for a given hVl and iVmax into the indi- 
vidual IR and UV contributions. This is done in Fig. 17 
for the ^Li fit of Fig 



15 



This figure verifies our prior 
claim that IR and UV corrections are roughly equal at the 
variational minima, while the IR(UV) correction rapidly 
dominates when we move to the right (left) of a minimum. 
If there are enough (/iJl, TVmax) points calculated where 
one of the two corrections is numerically insignificant, a 
simpler extrapolation with only three fit parameters is 
possible. An example of an IR-only fit is given for the 
NN-only calculation of ^"B at A = 2fm"^ IS given m 
Fig. [TSj Such IR-only fits are not possible for the cal- 
culations here including NNN because we do not have 
enough points sufficiently removed from the minimum in 

nn. 

The points in Fig. [18] were chosen to the right of the 
variational minimum (see Fig. [s]), where a fit to Eq. ^ 
implies that the UV correction to these points is much 
smaller than the IR correction. The observation that the 
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FIG. 18. (color online) Ground-state energy of ^"B for 
the NN-only evolved Hamiltonians at A = 2.0 fm~^ for 
iV^^„ = 4-8 and hO, = 28-32 MeV plus iVmax = 10 and 
hfl = 26, 28 MeV with an IR-only fit using Eq. Q. The fit 
value of i?oo is the dashed line. 



points with different hfl and V^ax values all lie on the 
same curve verifies that L is the correct variable. We 
note similar demonstrations for ^He results in Fig. 4 of 
Ref. [34] with the same NN-only interaction and in Fig. 11 
of Ref. [3F using a different NN interaction. The fit to 



E{L) ^Eoo+ B2e 



(6) 



for the energies shown in Fig. 18 is very good. However, 
the prediction for Eoo is about 0.6 to 1 MeV more bound 
than that from the "Extrapolation B" analysis. 

It would be premature to draw robust conclusions on 
the relative efficacy of the extrapolation schemes used 
here. In particular, further comparisons are needed 
where large A^max results are available to check small 
-Vmax extrapolations. However, for our present purposes 
it is sufficient that the results of the different schemes 
are consistent with each other to within the assessed 
uncertainties currently available. For the remainder of 
this work, we will use extrapolation procedures based on 
Eq. ([3|. 



IV. EVOLUTION 



A. Running of ground-state energies 



In Figs. [T9}|22| we show the dependence on A of the 
ground-state energy in '^Li, ^Be ^°B, and ^^C using our 
best estimate for infinite-basis space results based on 
"Extrapolation B" described in the previous section. For 
A < 1.5 fm~^ there is good convergence with small error 
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FIG. 19. (color online) Extrapolated ground-state energy of 
^Li as a function of A for each SRG calculation. The initial 
interaction was N^LO NN [2T] included up to iVA2max = 300 
and N^LO NNN [H |^ up to A^ASmax = 40. The dashed 
curves connect data points and error bars obtained using the 
extrapolations described in the text. The small black arrow 
on the left shows the experimental value. 
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FIG. 21. (color online) Same as Fig. 19 but for the B ground 



state. The small black arrow on the left shows the experimen- 
tal value. 
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FIG. 22. (color online) Same as Fig. 19 but for the C ground 
state, as well as three times the ^He ground state energy (solid 
curves). 



FIG. 20. (color online) Same as Fig. 19 but for the Be ground 
state, as well as twice the ''He ground state energy (solid 
curves). 



bars, although the numerical accuracy of the extrapo- 
lation with 3NFs is limited by the iVASmax cutoff. For 
A > 1.5 fm~^ the two criteria for selecting the optimal 
Ml lead to (slightly) different results, and their differ- 
ence grows with A; the shaded symbols in Figs. [19]- 
|22| correspond to using the Ml value that minimizes 
E{N-caayi) — £'00, whercas the open symbols correspond 



to using the Ml value that minimizes the numerical error 
estimate. For many cases we currently do not have data 
available to perform extrapolations at or near the Ml 
value that would minimize the numerical error estimate. 

Nevertheless, we see in these figures that for A between 
1.0 and 2.0 fm"^ the general pattern is the same as that 
observed for ''He in Fig. [l] We expect that the NN-only 
(green circles) start for A = oo at an (underbound) en- 
ergy but the larger A's required to verify this are not suf- 
ficiently converged here. However, the characteristic dip 
due to omitted induced NNN forces is clear in each of the 
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nuclei. Including the induced NNN matrix elements (blue 
squares) significantly reduces but does not eliminate the 
dependence on A. (Note that, as observed in Fig. [l] the 
NN-only curve should be equal to the NN+NNN induced 
result at large A, which is plausible from Figs. [T9}|22l but 
not verifiable in the present calculations.) The trend of 
the results with induced NNN interactions is increased 
binding as A decreases from 1.8 to 1.0 fm~^, consistent 
with ''He shown in Fig.[l] The magnitude of the decrease 
is about 0.7-1.1 MeV, without a systematic dependence 
on the nucleus. When initial NNN interactions are in- 
cluded (red diamonds), the qualitative dependence on A 
over this same range is similar, but the magnitude of the 
decrease is systematically larger by roughly 1 MeV or a 
factor of two (less in ^°B). The downward shift in ground- 
state energies takes them below the experimental values 
for the four nuclei of Figs. 19 22 The additional bind- 



ing provided by the initial 3NF increases from less than 
one MeV per nucleon in ''Li to more than 1.5 MeV per 
nucleon in ^^C, almost independent of A for A between 
1.0 and 2.0fm"\ 



The shape of an evolution curve (as in Figs. 19 -22) 
is determined by the interplay of short- and long-range 
effects in a given A-body sector. This has been demon- 
strated explicitly in Ref. [28] for a model but needs to 
be more systematically verified in realistic systems. For 
a given Hamiltonian (i.e. fixed initial interactions and 
fixed truncation of the evolution equations at the 2-body 
or at the 3-body level) we observe close similarity be- 
tween the various nuclear ground-state evolution curves 
presented here. Focusing on the range in A where the 
error bars are smaller (A = 1.0-1.8 fm~^), the energy 
variations for the NN-only calculations are reduced by a 
factor of 4-5 for the NN-I-NNN calculations (and more for 
NN-hNNN- induced). Thus the induced NNN interaction 
acts to (almost) restore A-independence in this region, 
with the residual variation attributed to four-body (and 
higher) forces. 

This suggests that the induced NNN interaction, as in 
the ^He case, is the leading correction to the NN-only 
results arising from SRG evolution. Moreover, the sim- 
ilar shapes of the evolution curves through the range of 
A is also consistent with the preservation of hierarchical 
induced many-body forces; that is, we expect that the 
induced four-body forces will provide a residual contri- 
bution smaller than the induced NNN interaction. This 
conclusion is also consistent with previous analysis of ex- 
pectation values for components of the evolved interac- 
tion [IS]. 

The net induced 4NF contribution in ^He was found 
to be a few hundred keV at A 1 fm^^ (see Fig. [I]). In 



Figs. 19 22 the effect of omitted induced four-body (and 
higher) forces at lower A (judging solely from the verti- 
cal range of each curve) is at least of order 1 to 2 MeV, 
depending on whether initial 3NFs are included. This is 
still small enough that for light nuclei it may be possible 
to exploit very highly evolved Hamiltonians, especially 
if a simple approximation for the 4NF contribution can 



be found. Alternatively, incorporating the induced 4NF 
contribution explicitly may become computationally fea- 
sible in the near future. A quantitative understanding 
of the magnitude and scaling of (assumed) induced 4NFs 
with A and with A is still lacking. This has become an 
important issue in light of the growing overbinding with 
larger A observed by Roth et al. when an initial 3NF is 
included (unless the 3NF cutoff is significantly reduced 
with respect to the NN cutoff) [T5] . 

In Fig. [20j we also show the evolution curves for two 
a-particles, noting that the A-dependence is significantly 
stronger for ^Be than for two a-particles. With the NN- 
only interaction, ®Be is actually bound for A < 2fm~^, 
but once the induced 3NF is included, it appears to be 
unbound by about 1 MeV. With the initial 3NFs wc find 
that ^Be is bound for A < 1.5 fm~^, but for larger values 
of A we cannot draw a firm conclusion. Experimentally, 
*Be is unbound by about 0.1 MeV (i.e. the lowest lying 
narrow resonance is 0.1 MeV above the two-a threshold). 
In Fig. [22] we also show three times the ground state en- 
ergy of^Ie. With the NN-only potential, "'^C is bound 
by about 20 to 25 MeV relative to three a particles with 
the same interaction for A between 1 and 2 fm~^, but once 
the induced 3NFs are taken into account, it is bound by 
only about 4 MeV. Once the initial N2L0 chiral 3NFs are 
incorporated, the binding relative to three a's increases 
to about 14 MeV, and is nearly independent of A. Exper- 
imentally, this energy difference is about 7.5 MeV. We 
expect that the overbinding by almost a factor of two 
relative to the three a threshold with the chiral N2L0 
3NFs will have important consequences for the low-lying 
spectrum of ^^C, in particular for the Hoyle state. 



B. Low-lying excited states 



In Figs. 23 - 25 we show the excitation energies of the 
lowest excited states in ^Li, ^Be, and ^"^B as function of 
hfl for different values of A at fixed A^max — 8. For ^Li 
and ®Be, the excitation energy of the first excited state is 
almost independent of both A and Ml^ and in good agree- 
ment with experiment. This independence suggests that 
these excitation energies are close to being converged. 
The second excited state of these nuclei shows a slight 
dependence on both A and HQ. In particular, at larger 
values of A the excitation energies show a variation with 
hfl indicating that these excitation energies are not yet 
converged. However, this variation is significantly less 
than the estimated extrapolation error in the absolute 
(ground state) energies for these nuclei. More interesting 
is that the (albeit small) dependence on A is significantly 
larger with initial 3NFs than with induced 3NFs only, as 
is evident by the smaller spread of the curves in the left 
panels of Figs. [23] and [24] 

In ^°B the situation is much more complicated. In 
Fig. [25] we show the excitation energies of the lowest 
four excited states of ^"B relative to the (3"'",0) (which 
is the ground state of ^"B): two (1^,0) states, a (2+,0) 
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state, and a (O"*", 1) state (analog state of the ground state 
of ^°Be and ^°C). It is immediately obvious that there is 
a much larger spread in the excitation energies than for 
^Li and ^Be, in particular with initial 3NFs. 

With induced 3NFs only, the lowest state is actually 
a (1~'',0) state, rather than the (3"'",0). This state how- 
ever is rather poorly converged relatively to the (3"'",0) 
state, as is evident from the strong Ml dependence of its 
(mostly negative) excitation energy, even for very small 
values of A. Furthermore, the second (1"'",0), as well as 
the lowest (2+, 0) and (0+, 1) states are very close to each 
other, all with excitation energies between 1 and 2 MeV, 
and all with a similar weak A and moderate Ml depen- 
dence. Experimentally, these three states have excitation 
energies of 2.154 MeV, 1.740 MeV, and 3.587 MeV respec- 
tively. Thus, in the absence of initial 3NFs, the chiral 
interactions not only predict the wrong ground state for 
but also a much too dense spectrum for the other 
low-lying states. 

With initial 3NFs the spectrum looks quite different. 



One of the two (1^,0) becomes strongly dependent on 
A, whereas the other (1^,0) state remains almost inde- 
pendent of A. In the region of parameter space where 
these two states are well-separated, they can be distin- 
guished by their quadrupole and magnetic moments. Al- 
though the quadrupole moments are not converged, they 
are clearly different for these two states: one has a small 
positive quadrupole moment of the order of one e fm^ or 



smaller (open symbols in Fig. 25), whereas the other has 
a negative quadrupole moment around —2 to —3 e fm^ 
(solid symbols in Fig. 25 ). The latter of these two states 
depends only weakly on A and Ml, and appears to be 
reasonably well converged, whereas the former is very 
strongly dependent on both A and Ml, and is not con- 
verged at all. However, for A from 1.5 fm^^ to 2.0 fm^^ 
these two states show significant mixing. (Experimen- 
tally, they are separated by about 1.4 MeV.) For com- 
parison, the quadrupole moment of the ground state is 
about +1 to -1-8 e fm^ with induced 3NFs only, and about 
-|-6 to +7 e fm^ with initial 3NFs, in reasonable agree- 
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ment with the experimental value of +8.472(56) e fm^, 
given the fact that the quadrupole moments are not yet 
converged in these basis spaces. It vifill take a major effort 
to develop robust extrapolation and error quantification 
tools for long-range operators, such as quadrupole opera- 
tors. We therefore defer a detailed consideration of these 
operators until those tools are developed. 

For small values of A, the slowly converging (!"*", 0) 
state with a small quadrupole moment actually becomes 
the lowest state, even with initial 3NFs. It is likely that 
for this state the induced four-body (and higher-body) 
interactions, which have been omitted in the present cal- 
culations, are important, though without convergence it 
is hard to draw firm conclusions. It is also possible that 
the strong A dependence is caused by the lack of conver- 
gence, and that once convergence (i.e. independence of 
both Njnax and hfl) is reached the results will be much 
less dependent on A. 

The other excited states shown in Fig. [25j namely the 
(0+,l) (open symbols) and the (2+,0) (solid symbols), 
show a strong dependence on both A and HQ with initial 
3NFs. Clearly, these excitation energies are not very well 
converged, but nevertheless we can see that the overall 
effect of the initial 3NFs is to increase their excitation 
energy and to separate these two states from each other. 



in qualitative agreement with the data. 

A striking difference between Figs. [23] and [24] on the 
one hand, and Fig. [25] on the other, is the strong A and 
hn dependence of the excitation energies in Fig. 25 (with 
the possible exception of one of the two low-lying(l+, 0) 
states) compared to the independence of A and hH. of the 
excitation energies in Figs. [23] and [24] One possible ex- 
planation for this observation is that the excited states 



in Figs. 23 and [24] can be interpreted as rotational exci- 
tations of the ground state jlD]. Thus these states have 
a very similar structure, and are therefore likely to ex- 
hibit a similar convergence pattern (i.e. dependence) 
and A dependence. We also see this in -'^^C, where the 
lowest 2+ and 4+ states form a rotational band with the 
ground state [33], but other low- lying states in ^^C are 
much more sensitive to the 3NFs and have a different 
convergence pattern 37J. 

Finally, we also calculated the magnetic moments of 
^Li, ''Be, and ^"B, see Table III The quoted numerical 
uncertainty in Table |III| includes both the dependence 
on the basis space parameters (i.e. hfl and iVmax depen- 
dence) and the A dependence. At iVmax = 8 the magnetic 
moments are typically converged to within a few percent, 
and the influence of the SRG evolution is less than a 
few percent, except for the | state of ^Be and the two 
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TABLE III. Magnetic moments for ''Li, ''Be, and with 
the NN+NNN-induced and NN+NNN interactions. 



Nucleus 


state 


induced 


IN iM iM 


expt . 


'^T i 


3 - 

2 








''T i 


1 - 
2 


— U.»(lj 






IjI 


7- 
2 


0.0[Z ) 


o.Zyo ) 




'Be 


3 - 
2 


-1.15(5) 


-1.15(6) 


-1.399 


'Be 


1 - 
2 


1.18(3) 


1.22(3) 




'Be 


7 - 
2 


0.24-0.56 


0.37-1.07 




lOg 


3+ 


1.85(1) 


1.83(2) 


1.8006 


lOg 


1+ 


0.84(2) 


0.78-0.85 


0.63(12) 


lOg 


1+ 


0.35(2) 


0.34-0.41 





(l+,0) states in -'^^B. A closer look at the different com- 
ponents contributing to the magnetic moment of the | 
state of ^Be shows that the contributions from the neu- 
tron intrinsic spin and from the proton angular momen- 
tum nearly cancel, leaving the proton intrinsic spin con- 
tribution to dominate the magnetic moment, both with 
and without initial 3NFs. For the two states in ^°B, most 
of the parameter dependence of the magnetic moments 
for these two states is due to (strong) mixing. In general, 
adding initial 3NFs to the chiral N3L0 NN interaction 
does not have a significant effect on these magnetic mo- 
ments. 

These magnetic moments are calculated in impulse 
approximation, using the canonical Ml operator. Of 
course, we should use a current operator that is consis- 
tent with the chiral Hamiltonian that we are using: for 
the ground state of ^Li and ^Be we might expect, based 
on Refs [HI |35] , a correction due to meson exchange cur- 
rents of about 10%, in the direction that would bring our 
results toward agreement with experiment. 



V. CONCLUSIONS 

We have presented ab initio NCFC calculations of en- 
ergies in the p-shell using SRG-evolved two- and three- 
nucleon forces. Several different procedures were con- 
sidered to extrapolate energies to infinite harmonic- 
oscillator basis size; for the range in the evolution pa- 
rameter A we focus on (from 1 to 2fm~^) they give 
consistent results within estimated error bars. Error 
bars above 2fm~^ grow rapidly and limit what we can 
quantitatively conclude about A dependence. As antic- 
ipated from results with lighter nuclei, inclusion of in- 
duced NNN interaction significantly reduced the A de- 
pendence of ground-state energies compared to NN-only 
calculations for all nuclei considered. Furthermore, the 
rapid improvement of convergence with decreasing A of 
ground-state and low-lying state energies observed for 
NN-only calculations |14j carries over when only induced 
NNN interactions are included. With initial NNN inter- 
actions, ground-state convergence is similarly improved 



but some excited states show very different behavior, see 
Fig. [251 

With NN+NNN-induced interactions (but without ini- 
tial NNN interactions), the net change in total ground- 
state energy for nuclei in the p-shell was found to be 
small (and within extrapolation error bars) between A = 
2.0 fm~^ and A = 1.5 fm""'^, but systematically decreases 
(becomes more bound) as A decreases to 1.0 fm^^ by 
about 1 MeV. The A dependence is small. This addi- 
tional binding can be attributed to four- and higher-body 
forces, which is of natural size (as implied by EFT power 
counting) despite the extreme degree of softening. While 
this might appear to be a narrow range in A, we empha- 
size that there is significant evolution (e.g., note that the 
natural SRG evolution variable s = l/A** increases by a 
factor of 16 as A decreases from 2 to 1 fm^^). When initial 
3NFs are included, a similar pattern is found, but the de- 
crease is more nucleus dependent and larger by as much 
as a factor of 2 (i.e., up to about 2 MeV). This increase 
in binding is consistent with the difference in running 
between NN-hNNN-induced and NN-^NNN observed for 
^^C in Ref. [18 with the same interactions. 

Results with 3NFs for higher A are consistent with 
small changes, but the uncertainties after extrapolation 
are too large to be definitive. In contrast. Roth et al. 
found a steady linear increase in binding from A above 
2.2 fm"^ down to 1.6fm"^ for [18|. This is not evi- 
dent in the systematics of the central values of our extrap- 
olations, but possible within the estimated error bars. It 
would be most helpful to have accurate energy calcula- 
tions for the initial Hamiltonian to fully assess the degree 
of running down to A = 2 fm~^ . For excited states we also 
find similar quantitative differences between NN-I-NNN- 
induced and NN+NNN calculations and for some partic- 
ular states there are qualitative differences. 

Further investigations are warranted. For example, it 
will be important to compare our present results, which 
use harmonic oscillator evolved SRG 3NFs, to forth- 
coming calculations using the same initial interactions 
but evolved in momentum space using recently devel- 
oped SRG technology [20]. Further explorations at low 
A will help to map out the quantitative scaling of in- 
duced 4NF contributions. Improved convergence at these 
low resolutions motivates searching for simple approxi- 
mations to account for these 4NF contributions to the 
energy and other observables. Besides additional CI cal- 
culations, application of highly evolved two- and three- 
nucleon forces may also be fruitful for coupled cluster [JS] 
methods, ab initio density functional theory [441 , and 
NCSM/RGM 5, 45 48 calculations of light nuclear re- 
actions. 

Our studies were limited in basis size by the avail- 
able codes and computer resources. We anticipate fur- 
ther developments in improved basis construction and 
evolution algorithms. We could also study non-oscillator 
basis spaces [49l [5^ and apply the importance trunca- 
tion technique [T7| to increase available basis sizes. We 
look forward to a detailed extrapolation analysis of such 
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results. In addition, a difficulty encountered here with 
an insufficient initial A = 3 basis at smaller Ml (see Sec- 
tion III A) might be circumvented by the momentum- 



space evolution technology [5D], so that matrix elements 
in the oscillator basis are only calculated after evolution. 

We considered only a single initial NN-fNNN Hamil- 
tonian in our analysis. While in past investigations jl6j 
other choices of initial NN Hamiltonians have displayed 
no qualitative difference in the effects of the SRG pro- 
cedure, studies with a range of Hamiltonians are highly 
desirable. It will also be important to consider consis- 
tent operators for other observables. In this regard, new 
extrapolation methods for radii and other long-distance 
operators may be particularly valuable [34j . 

Our results demonstrate a level of precision in induced 
many-body force effects that allows further analysis of 
the impact of additional xEFT inputs. In all cases pre- 
sented here the addition of 3NFs from the initial xEFT 
Hamiltonian overbinds the ground states. Our results 
suggest that the effects of missing induced forces due to 
softening transformations are small (at least in larger nu- 
clei) compared to discrepancies with experiment. These 
discrepancies may be reduced by additional 3NF and 4NF 
contributions at N^LO in the xEFT. Suitable matrix el- 
ements of these terms will be available for calculations in 
the near future and it will be important to incorporate 



them. 



ACKNOWLEDGMENTS 

This work was supported in part by the National Sci- 
ence Foundation under Grant Nos. PIIY-1002478 and 
PHY-0904782 and the Department of Energy under 
Grant Nos. DE-FG02-87ER40371, DE-FC02-07ER41457 
(SciDAC-2/UNEDF), DE-FC02-09ER41582 (SciDAC- 
2/UNEDF) and DESC0008485 (SciDAC-3/NUCLEI). 
Prepared in part by LLNL under Contract DE-AC52- 
07NA27344. Support from the Natural Sciences and En- 
gineering Research Council of Canada (NSERC) Grant 
No. 401945-2011 is acknowledged. TRIUMF receives 
funding via a contribution through the National Research 
Council Canada. A portion of the computational re- 
sources were provided by the National Energy Research 
Scientific Computing Center (NERSC), which is sup- 
ported by the DOE Office of Science, and by an INCITE 
award, "Nuclear Structure and Nuclear Reactions", from 
the DOE Office of Advanced Scientific Computing. This 
research also used resources of the Oak Ridge Leader- 
ship Computing Facility at ORNL, which is supported 
by the DOE Office of Science under Contract DE-AC05- 
00OR22725. 



[1] P. Maris, J. P. Vary, and A. M. Shirokov, Phys. Rev. C 
r 79, 014308 (20 09) , arXiv:0808.3420 [nucl-th]^ 
[2J P. Navratil, S. Quaglioni, I. Stetcu, and B. R. Barrett, 
'J. Phys. G 36, 083101 (2009) , arXiv:0904.0463 [nucl-th]l 
]3] P. Maris, A. Shirokov, and J. Vary, Phys. Rev. C 81^ 
021301 (2010), arXiv:0911.2281 [nucl-th]. 



and D. Dean, 


Phys. Rev. Lett. 106, 202502 (2011) 


,arXiv:1101.5124 [nucl-th| 





5]^?. Navratil and S. Quaglioni, | Phys. Rev. Lett. 108 
042503 (2012), arXiv:1110.0460 [nucl-th] , 
B. R. Barrett, P. Navratil, and J. P. Vary, Prog. Part 
"Nucl. Phys. 69, 131 (2013) 



[7] S. D. G lazek and K. G. Wilson, ,Phys. Rev. D 48, 5863 
(1993)1 



[8] F. Wegner, Ann. Phys. 506, 77 (1994). 
[9] S. K. Bogner, R. J. Furnstahl, and A. Schwenk, Prog. 
Part. Nucl. Phys. 65794' (20l0)| |arXiv:0912.3688 [nucl- 



ih]r 



[10] R. Furnstahl, Nucl. Phys. Pr oc. Suppl. 228, 139 (2012) 

[arXiv: 1203.1779 [nucl-th] 
[11] S. K. Bogner, R. J. Furnstahl, and R. J. Perry, Phys. 

Rev. C 75, 061001 (2007), nucl-th/0611045 
[12] S. K. Bogner, R. J. Furnstah l, R. J. Perry, and 

A. Schwenk, |Phys. Lett. B 649, 488 (2007)[ |arXiv:nucl-| 

[ ~th/07 01013"', 

[13] E. D. Jur genson, S. K. Bogner, R. J. Furnstahl, and R. J 



Perry, Phy s. Rev. C 78, 014003 (2008)| |arXiv:0711.4252 

[nucl-th] 

14] S. K. Bogner, R. J. Furnstahl, P. Maris, R. J. Perry, 
A. Schwenk, and J. P. Vary, Nucl. Phys. A 801, 21, 



[18 
[19 

[2o; 

[21 
[22; 



[23 
[24 



(2008) 



I arXiv:07 08.3754' [n ucl-th]] 



E. D. Jurgenson, P. Navratil, and R. J. Furnstahl, 'PhysJ 
Rev. Lett. 103, 082501 (2009) , arXiv:0905. 1873 [nucl-tl^ 
E. Jurgenson, P. Navratil, and R. Furnstahl, Phys. Rev!] 
"C 83, 034301 (2011), arXiv:1011.4085 [nucl-th], 
R. Roth, Phys. Rev. C 79, 064324 (2009)| 
arXiv:0 903.4605 [nucl-th], 

R. Roth, J. Langhammer, A. Calci, S. Binder, and 
P. Nav ratil, 'Phys. Rev. Lett. 107, 072501 (2011)| 
|arXiv:1105.3173 [nucl-th] 



n. rtorn, d. oina 
and P. Navratil, 


:;r, IX. voDig, A. i^aici, j. ijangnammer 
Phys. Rev. Lett. 109, 052501 (2012) 


|arXiv:1112.0287 


nucl-th] 



Rev. C 85, 021002 



(2012)1 



K. Hebeler, 
|arXiv: 1 20 1.0169' [nucl-th] 

D. R. Ente'm and R. Machleidt, Phys. Rev. C 68, 041001 
(2003), nucl-th/0304018 

E. Epelbaum, A. Nogga, W. Gloeckle, H. Kamada, U.- 
G. Meissner, an d H. Witala, Phys. Rev. C 66, 06400T| 
(2002)1 |arXiv:nu cl-th/0208023 , 
P. Navratil , ^Few Body 
|arXiv:0707.4680 [nucl-th] ~ 

P. Navratil, V. G. Gueorguiev, J. P. Vary, W. E. Ormand, 
and A. Nogga, Phys. Rev. Lett. 99, 042501 (2007), nucl- 
th/0701038 

Kehrein, The Flow Equation Approach to Many- 
Particle Systems (Springer, Berlin, 2006). 
E. Anderson et ai, Phys . Rev. C 77, 037001 (2008)| 
arXiv:0801.1098 [nucl-th] 



Syst. 41, 117 (200f)] 



W. Li, E. Anderson, and R. Furnstahl, Phys. Rev. C 84, 
054002 (2011), arXiv:1106. 2835 [nucl-th]. 



19 



[28] E. D. Jurgenson and R. J. Furnstahl, N ucl. Phys. A 818, 
r 1 52 (2009) , arXiv:0809.4199 [nucl-thj^ 



[29j D . Gazit, S. Quaglioni, and P. Navratil, Phys. Rev. Lett 
f 103, 102502 (2009), arXiv:0812.4444 [nucl-thJJ 
[30] P. Sternberg, E. G. Ng, C. Yang, P. Maris, J. P. Vary, 
M. Sosonkina, and H. V. Le, in Proceedings of th3{ 
~2008 ACM/IEEE conference on SupercompuUng , SC '08 



(IEEE Press, Piscataway, NJ, USA, 2008) pp. 15:1-15:12. 
[31] P. Maris, M. Sosonkina, J. P. Vary, E. G. Ng, and 

C. Yang, Procedia CS 1, 97 (2010). 

[32] H. M. Aktulga, C. Yang, E. G. Ng, P. Maris, and 
J. P. Vary, in Euro-Par, Lecture Notes in Computer Sci- 
ence, Vol. 7484, edited by C. Kaklamanis, T. S. Pap- 
atheodorou, and P. G. Spirakis (Springer, 2012) pp. 830- 
842. 

[33] P. Maris, H. Aktulga, M. Caprio, U. Qatalyiirek, E. Ng, 

D. Oryspayev, H. Potter, E. Saule, M. Sosonkina, J. Vary, 
C. Yan g, and Z. Zhou, |J. Phys. Conf. Ser. 403, 0120191 
(2012)1 



[34] R. Furnstahl, G. Hag 


en, and T. Papenbrock, Phys. Rev. 


C 86, 031301 (2012) 


arXiv: 1207.6100 [nucl-th] 





P. Maris, and J. P. Vary,' Phys. Rev. C 86, 054002 (2012), 

arXiv:1205.3230 [nucl-th]', 
[36] S. More, A. Ekstrom, R. J. Furnstahl, G. Hagen, and 

T. Papenbrock, (2013), arXiv : 1302.3815 [n ucl-th] 
[37] A. Calci, J. Langhammer, P. Maris, R. Roth, and 

J. Vary, (2013), in preparation. 



C. Cockrell, J. P. Vary, and P. Maris, P hys. Rev. C 86, 
_034325 (2012), ,arXiv:1201. 0724 [nucl-th]^ 

E. Epelbaum, H.-W. Hammer, and U.-G. Meissner, Rev.| 
"Mod. Phys. 81, 1773 (2009)', 'arXiv:0811.1338 [nucl-th]' 



M. C aprio, P. Maris, and J. Vary,jPhys. Lett. B719, 1791 



(2013), arXiv:1301.0956 [nucl-th] 
L. E. Marcucci, M. Pervin, S. C. Pieper, R. Schiavilla, 
and R. B. Wiringa, Phys. Rev. C 78, 065501 (2008) 
S. Pastore, S. C. Pieper, R. Schiavilla, and R. Wiringa, 
(2012), | arXiv:1212. 3375 [nucl-t h] 

G. Hagen, D. J. Dean, M. Hjorth- Jensen, T. Papen- 
brock, and A. Schwenk, Phys. Rev. C 76, 044305 (2007)^, 
arXiv:0707.1516 [nucl-th] 
J 



E. 



Prut and L. Platt er, |Phys. Rev . C 84, 014318 
(2011)1 |arXiv:110 4.4357 [n ucl-thj | 

^ Oiiacrlinni nnri P IVaArra + il H 



Lett. 101, 



S. Quaglioni and P. Navratil, 'Phys. Rev 
"092501 (2008), arXiv:0804.1560 [nucl-th] | 
P. Navratil, R. Roth, and S. Quaglioni, Phys. Rev. C 
82, 034609 (2010), arXiv: 1007.0525 [nucl-thj^j^ 



P. Navratil, R. Roth, and S. Quaglioni, P hys.Lett. B704, 
379 (2011), arXiv:1105.5977 [n ucl-th] 
S. Quaglioni, P. Navratil, R. Roth, and W. Horiuchi, 



(2012), arXiv:1203.02 68 [nucl-th]| 

M. Caprio, P. Maris, and J. Vary, |Phys. Rev. C 86, 



034312 (2012), arXiv:1208.4156 [nucl-thf 
M. Caprio, P. Maris, and J. Vary, J. Phys. Conf. Ser. 



403, 012 014 (2012)1 |arXiv: 1208^6055 [nucl-thJ 



